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ABSTRACT 

In this work, we consider the cosmological constraints on the interacting dark energy models. We 
' generalize the models considered previously by Guo et al. p^ . Costa and Alcaniz '16], and try to 

^SJ \ discuss two general types of models: type I models are characterized by p^j pm = f{a) and fia) can 

be any function of scale factor a, whereas type II models are characterized by p m — PmO ^ 
and e{a) can be any function of a. We obtain the cosmological constraints on the type I and II 
models with power-law, CPL-like, logarithmic f{a) and e(a) by using the latest observational data. 
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I. INTRODUCTION 



The dark energy has been one of the most active fields in modern cosmology since the discovery of the 
accelerated expansion of our universe (see e.g. [l| for reviews) . Among the conundrums in the dark energy 
cosmology, the so-called cosmological coincidence problem is the most familiar one. This problem is asking 
why are we living in an epoch in which the densities of dark energy and matter are comparable? Since 
their densities scale differently with the expansion of our universe, there should be some fine-tunings. To 
alleviate the cosmological coincidence problem, it is natural to consider the possible interaction between 
dark energy and dark matter in the literature (see e.g. [2-9, 31, 32]). In fact, since the nature of both 
dark energy and dark matter are still unknown, there is no physical argument to exclude the possible 
interaction between them. On the contrary, some observational evidences of this interaction have been 
found recently. For example, in a series of papers by Bertolami et al. [Toj . they shown that the Abell 
Cluster A586 exhibits evidence of the interaction between dark energy and dark matter, and they argued 
that this interaction might imply a violation of the equivalence principle. On the other hand, in 
Abdalla et al. found the signature of interaction between dark energy and dark matter by using optical. 
X-ray and weak lensing data from 33 relaxed galaxy clusters. Therefore, it is reasonable to consider the 
interaction between dark energy and dark matter in cosmology. 

We consider a flat Friedmann-Robertson- Walker (FRW) universe. In the literature, it is usual to assume 
that dark energy and dark matter interact through a coupling term Q, according to 

Pm + SHprn = Q, (1) 

p^+3Hp,il+w,)^ ^Q, (2) 

where pm and p^ are densities of dark matter and dark energy (we assume that the baryon component 
can be ignored); is the equation-of-state parameter (EoS) of dark energy and it is assumed to be a 
constant; a dot denotes the derivative with respect to cosmic time t\ H = a/a is the Hubble parameter; 
a = (1 -I- z)~^ is the scale factor (we have set oq = 1; the subscript "0" indicates the present value 
of corresponding quantity; z is the redshift). Notice that Eqs. (HJ and ^ preserve the total energy 
conservation equation 

ptot + 3i7/9tot(l + Woff) = 0, (3) 

where ptot = Px ~^Pm is the total energy; Wos is the total (effective) EoS. Since there is no natural guidance 
from fundamental physics on the coupling term Q, one can only discuss it to a phenomenal level. The 
most familiar coupling terms extensively considered in the literature are Q = aK,pm4>i Q = 'iPHptot, 
and Q = 3r]Hpm- The first one arises from, for instance, string theory or scalar-tensor theory (including 
Brans-Dicke theory) [1-[6| . The other two are phenomenally proposed to alleviate the coincidence problem 
in the other dark energy models 043- In the usual approach, one should priorly write down the coupling 
term Q, and then obtain the evolutions of pm and p^ from Eqs. ^ and ([2]), respectively. In fact, this is 
the common way to study the interacting dark energy models in the literature. 

However, there is also an alternative way in the literature [l^ - [l6j . One can reverse the logic mentioned 
above. Due to the interaction Q, the evolutions of pm and p^ should deviate from the ones without 
interaction, i.e., pm oc a"^ and p^ oc a~'^^^'^^x \ respectively. If the deviated evolutions of pm and/or 
Px are given, one can find the corresponding interaction Q from Eqs. ([1]) and ([2]). Naively, the simplest 
example has been considered by Wang and Meng p,2], namely 

Pm = Pmo a"^''"^ (4) 

where e is a constant which measures the deviation from the normal pm oc . Substituting into Eq. ([T]), 
it is easy to find the corresponding interaction Q — eHpm 0, [12, [13, [i3| ■ Alternatively, one can consider 
another type of interacting dark energy model which is characterized by [3, [Tsl 



— u , 

Pm PmO 



(5) 



where ^ is a constant which measures the severity of the coincidence problem. From Eqs. ((TJ, ^ and ([5]), 
one can find that the corresponding interaction is given by [l5j | 

Q - -Hp^ilx (e + ) = -Hp^ r!,„ (e + ^w^), (6) 
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where Hi = 8TTGpi/{3H^) for i — m and X, which are the fractional energy densities of dark matter 
and dark energy, respectively. In fact, Guo et al. [l5| considered the cosmological constraints on the 
interacting dark energy model characterized by Eq. ([5]) with the 71 SNLS Type la supernovae (SNIa) 
dataset, the shift parameter R from the Wilkinson Microwave Anisotropy Probe 3-year (WMAP3) data, 
and the distance parameter A of the measurement of the BAO peak in the distribution of SDSS luminous 
red galaxies. On the other hand, the interacting dark energy model characterized by Eq. ^ has been 
extended in [I^. It is more realistic that e is a function of time. Costa and Alcaniz [l^ considered the 
interacting dark energy model characterized by 

P™-p™oa"'+^^'^\ (7) 

in which e(a) was chosen to be 

e(a)=eoa'\ (8) 

where eo and ei are constants. They obtained the constraints on this model by using the 307 Union SNIa 
dataset, the CMB constraint Vlmoh?' = 0.109 ± 0.006 from the Wilkinson Microwave Anisotropy Probe 
5-year (WMAP5) data, and the distance ratio from zbao — 0.35 to z^s = 1089 measured by SDSS, 
namely Rbao/ls = 0.0979 ± 0.0036. 

In the present work, we generalize the interacting dark energy models considered in [l6j |. and we 
call them type I and II models, respectively. The type I models are characterized by 

^-/(a), (9) 

Pm 

where /(a) can be any function of a, beyond the special case in Eq. ([5]). The type II models are character- 
ized by Eq. ([7]) but e(a) can be any function of a, beyond the special case in Eq. (|5]). In the present work, 
we consider the constraints on these models by using the latest cosmological observations, namely, the 
397 Constitution SNIa dataset [l^, the sh ift p arameter R from the newly released Wilkinson Microwave 
Anisotropy Probe 7-year (WMAP7) data [I^, and the distance parameter A of the measurement of the 
BAO peak in the distribution of SDSS luminous red galaxies [l^, [2l|. In the next section, we briefly 
introduce these observational data. In Sec. |IlT] and Sec. llVi we discuss the type I and II models, and 
consider their cosmological constraints, respectively. A brief summary is given in Sec. |Vl 



II. OBSERVATIONAL DATA 



In the present work, we will consider the latest cosmological observations, namely, the 397 Constitution 
SNIa dataset [l^ , the shift parameter R from the newly released Wilkinson Microwave Anisotropy Probe 
7-year (WMAP7) data [l^, and the distanceparameter A of the measurement of the BAO peak in the 
distribution of SDSS luminous red galaxies 0, . 

The data points of the 397 Constitution SNIa compiled in [l^ are given in terms of the distance 
modulus Pobs{zi). On the other hand, the theoretical distance modulus is defined as 

Pth{zi) = b\ogif^DL[zi) + po 1 (10) 
where /zq = 42.38 — 51og]^Q h and h is the Hubble constant Hq in imits of 100 km/s/Mpc, whereas 

in which E = H/Hq, and p denotes the model parameters. Correspondingly, the from the 397 
Constitution SNIa is given by 
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where a is the corresponding \a error. The parameter /iq is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginahzation over /iQ. However, there is an alternative 
way. Following 22, 23], the minimization with respect to /io can be made by expanding the of Eq. (jl2p 
with respect to /io as 

X^(p) = i-2/.oB + Ai§C, (13) 

where 



X ^ f^obsjzt) - fithjzj; = 0,p) ^ ^ 1 



Eq. (I13p has a minimum for /iq = B/C at 



S(p) 



C 

^.2 



X^(p)=A(p)-^. (14) 



Since xfi^min — xfi^min obviously, we can instead minimize xfi which is independent of ^o- 

There are some other observational data, such as the observations of cosmic microwave background 
(CMB) anisotropy [l§] and large-scale structure (LSS) 20]. However, using the full data of CMB and LSS 
to perform a global fitting consumes a large amount of computation time and power. As an alternative, one 
can instead use the shift parameter R from the CMB, and the distance parameter A of the measurement 
of the baryon acoustic oscillation (BAO) peak in the distribution of SDSS luminous red galaxies. In the 
literature, the shift parameter R and the distance parameter A have been used extensively. It is argued 
that they are model-independent |24l] . while R and A contain the main information of the observations 
of CMB and BAO, respectively. 

As is well known, the shift parameter R of the CMB is defined by [U, [2^ 

where the redshift of recombination — 1091.3 which has been updated in the Wilkinson Microwave 
Anisotropy Probe 7- year (WMAP7) data [l^. The shift parameter R relates the angular diameter 
distance to the last scattering surface, the comoving size of the sound horizon at z* and the angular 
scale of the first acoustic peak in CMB power spectrum of temperature fluctuations [H, [2^ . The value 
of R has been updated to 1.725 ± 0.018 from the WMAP7 data On the other hand, the distance 
parameter A of the measurement of the BAO peak in the distribution of SDSS luminous red galaxies [l^l 
is given by 



1 f-" dS "1^/^ 



Zb Jo E{z) 



(16) 



where Zb = 0.35. In [U, the value of A has been determined to be 0.469 (ns/0.98)-°-35 ± 0.017. Here the 
scalar spectral index is taken to be 0.963, which has been updated from the WMAP7 data [l^. So, 
the total is given by 

^xl + XcMB + x\ao ' (17) 

where x^ is given in Eq. ([14]), Xcmb = {R - RobsY/<JR and x^bao = {A - AobsY/<j\- The best-fit 
model parameters are determined by minimizing the total x^- As in ^ 113], the 68% confidence level 
is determined by Ax^ = x^ — Xmin — 1-0' 2.3 and 3.53 for np = 1, 2 and 3, respectively, where Up 
is the number of free model parameters. Similarly, the 95% confidence level is determined by Ax^ = 
~ Xmin — 4-0, 6.17 and 8.02 for Up — 1, 2 and 3, respectively. 
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III. TYPE I MODELS 



A. Equations 

As mentioned in Sec. HI the type I models are characterized by Eq. ([9]), whereas /(a) can be any 
function of a. From Eq. it is easy to obtain 

= , a„ = . (18) 

Substituting = pmf{cL) into Eq. ([2]) and using p„j from Eq. ([T]), we can find that the corresponding 
interaction term is given by 

Q = -Hprn^^ ^fly + 3w^^ = -Hp^VL,n ^Oy + ^ , (19) 

where a prime denotes the derivative with respect to a. Obviously, if /(a) oc a^, Eq. (1191) reduces to 
Eq. © which has been obtained in [l^. On the other hand, one can recast Eq. ^ as 

^^^ = -3(l + Weff) = -3(l + l],u;,). (20) 
am a 

Using Eq. (|T8l) . we can integrate Eq. (l20l) to obtain 



Ptot ^ a ^ exp (^^ J J In ■ const. , (21) 

where const, is an integral constant, which can be determined by requiring the condition ptotia = 1) = 
Ptot,o = 3Hq / (SttG). Once ptot is on hand, we can readily find the corresponding E = H/Hq from 
Friedmann equation, and then fit it to the observational data. Correspondingly, p.^ ~ Q^ptot and 
Pm = ^mPtot are also available from Eqs. ^TE\\ and (l?T|) . Finally, it is worth noting that by definition 
we have 

/o = /(a ^ 1) = ^ = - 1 , (22) 

PmO i 'mO 

which is useful to fix one of parameters in /(a). 



B. Cosmological constraints on type I models 

In this subsection, we consider the cosmological constraints on type I models by using the observational 
data given in SecHT] At first, we consider the power-law case with 

/(a) = /oa«, (23) 

where ^ is a constant; /q can be determined by definition ^ to be the one given in Eq. ([2^ . and 
hence it is not an independent parameter. In this case, there are three free model parameters, namely, 
^mo, and ^. Although the cosmological constraints on the model characterized by Eq. has been 
considered by Guo et al. 15], as mentioned in Sec. HI they have used the earlier observational data. 
Therefore, it is still worthwhile to consider the cosmological constraints once more in the present work by 
using the latest observational data mentioned in Sec.|lll Substituting Eq. ([^51) into Eq. (PT|) and requiring 
Ptot{o- = 1) = Ptotfi, we can determine the integral constant, and finally obtain 

Ptot = Ptot,oa^^[^mO + il-^mo)a^] (24) 
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Substituting into Friedmann equation, we find that 



H 



(1 + zf [r!„o + (1 - r!™o) (1 + z)-^ ] '""^^ 



(25) 



By minimizing the corresponding total in Eq. (1171) . we find the best-fit parameters fimo — 0.281, 
= —0.982 and ^ = 2.988, whereas Xmin = 465.604. In Fig. [U we also present the corresponding 68% 
and 95% confidence level contours in w-^ — ^ plane, fimo — C plane and Q-mO — plane. It is easy to 
see that these constraints on the model characterized by /(a) = foa^ are much tighter than the ones 
obtained by Guo et al. [l^ , thanks to the latest observational data. 
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FIG. 1: The 68% and 95% confidence level contours in ui^ — C plane, fimo — £, plane and Q.mO — plane for the 
type I model characterized by /(a) = fo a^. The best-fit parameters are also indicated by the black solid points. 
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0.22 0.23 0.24 0.25 0.26 0.27 



FIG. 2: The same as in Fig. [T] except for the type I model characterized by /(a) — fo + ^{l — a) with the condition 
^ > 1 — f^mo- See the text for details. 

Next, we consider a new case with 

/(a) = /o+^(l-a), (26) 

which can be regarded as a Unear expansion of /(a) with respect to a, similar to the familiar Chevallier- 
Polarski-Linder (CPL) parameterization for EoS 'w{a) — wo + Wail — a) [28]. Again, /o can be determined 
by definition ^ to be the one given in Eq. (j22p . and hence it is not an independent parameter. Thus, 
there are three free model parameters, namely, ^Imo, w-^ and ^. Substituting Eq. (j26p into Eq. (j2ip and 
requiring ptot{o- = 1) = Ptot.Oi we can determine the integral constant, and finally obtain 

Ptot=Ptot,oa X ' [{1 + £,ilmoj a - t.ilm.0 \ ■ (27) 
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Substituting into Friedmann equation, we find that 

= (l + z)3(l+-x)[(l+eano)(l + ^)-Cano]"'"'-'''"''/''''^''""^ (28) 

Noting Eq. © and imposing the condition > 0, we have > 1 — il„Q from Eqs. (pS)) and (P^ . Under 
this condition, by minimizing the corresponding total in Eq. (jl7p . we find the best- fit parameters 
fi^o = 0.218, = -0.483 and ^ = -3.584, whereas x^,„ = 563.77. In Fig. d we also present the 
corresponding 68% and 95% confidence level contours in — plane, r2,„o — S, plane and il„jo — 
plane. It is easy to see that the 68% and 95% confidence level contours are very close. On the other 
hand, x^i„ = 563.77 is fairly larger than the degree of freedom dof ~ 400. 




-1 



0.28 0.285 0.29 0.295 

FIG. 3: The same as in Fig. [1] except for the type I model characterized by /(a) = /o + C(l ~ without the 
condition ^ > 1 — n~g. See the text for details. 
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So, we give up the condition ^ > 1 — f^mo case of f (a) = fo + £,{l — a). This means that might 

be negative in the early universe. In fact, Guo et al. [15] also explicitly include this possibility. Since 
such a negative energy appears in phantom models [29] and modified gravity models [30| , it is reasonable 
to consider this possibility. Without the condition ^ > 1 — ^^"q, by minimizing the corresponding total 
in Eq. (jl7p . we find the best-fit parameters — 0.288, w,^ ~ —0.868 and ^ — —2.957, whereas 
Xmin — 467.718. In Fig. [3l we also present the corresponding 68% and 95% confidence level contours in 
w-j^ — ^ plane, Q,m{) — S, plane and fimo ~ plane. Obviously, these results are significantly better than 
the ones with the condition ^ > 1 — ^l„o, whereas the corresponding x^i„ = 467.718 is also better. 

Finally, one might consider the logarithmic case with 

/(a) = ./o + ^lna, (29) 

which can be regarded as a linear expansion of /(a) with respect to the so-called e-folding time A/" = Ina 
in the literature. This case seems attractive since in Eq. (PT|) the integration is with respect to Ina. 
However, in this case, /(a) (and hence p^) diverges when a — > in the early universe. So, we do not 
consider the logarithmic case in type I model. 



IV. TYPE II MODELS 



A. Equations 

In this section, we turn to type II models, which are characterized by Eq. ([7]) whereas e(a) can be any 
function of a. Substituting Eq. ([7]) into Eq. (jlj, we can easily find that the corresponding interaction 
term is given by 

Q = Hp,n[e{a) + ae'{a)hYa]. (30) 

It is worth noting that in type II models, there is no condition like Eq. (1^^ in type I models to reduce 
the number of free model parameters. If there are at least two parameters in e(a), adding and w-^, 
we should have four free model parameters or even more. In this case, the constraints will be very loose, 
and the calculations will be very involved. Instead, we follow Costa and Alcaniz to consider the case 
with a fixed = —1, namely, the role of dark energy is played by a decaying A [H, [1^. Therefore, 
Eq. ([2]) becomes 

PK^-Q^-Hpm[t{a) + ae{a)\na]. (31) 
Substituting Eq. into Eq. ([3l|), we have 

^ = -p™o a-^+<'''> [ e(a) + ae' [a) In a ] . (32) 
da 

We can integrate Eq. to obtain 

PA=PmO a~*+''^°-^e{d) + ae'{a)lna]da + pAo ■ (33) 



Substituting Eqs. ([55)1 and ([7]) into Friedmann equation, we find that 

ij2 



where 



E^ = ^-^ = n,„o 0{a) + (1 - r!,„o) , (34) 

^0 



r-l 

(a) = a-3+^('') -f / &-^+<^^e{d) + &e'{d)lnd]d&. (35) 
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B. Cosmological constraints on type II models 



In this subsection, we consider the cosmological constraints on type II models by using the observational 
data given in Sec.|lll At first, we consider the power-law case with 

e(a)=eoa'\ (36) 

where eo and ei are constants. In this case, there are three free model parameters, namely, QmOi Sind ei. 
Although the cosmological constraints on the model characterized by Eq. (j36p has been considered by 
Costa and Alcaniz [l^, as mentioned in Sec. HI they have used the earlier observational data. Therefore, 
it is still worthwhile to consider the cosmological constraints once more in the present work by using the 
latest observational data mentioned in SecHH Substituting Eq. into Eqs. and (|55|) . we can then 
fit this model to the observational data. By minimizing the corresponding total in Eq. (|17p . we find 
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the best-fit parameters ilmo = 0.282, eo = —0.129 and ei = 1.263, whereas Xmm = 465.635. In Fig. |4l we 
also present the corresponding 68% and 95% confidence level contours in eo — ei plane, ilmo — eo plane 
and flmQ — ei plane. It is easy to see from the cq — ei plane that if eo is close to zero, ei cannot be too 
negative. On the other hand, the constraint on the parameter ei is still very loose. 
Next, we turn to the CPL-like case with 

e(a) = eo + ei(l-a), (37) 

which can be regarded as a linear expansion of e(a) with respect to a, similar to the well-known CPL 
parameterization for EoS w{a) = wq + Wa{l — a) [28j]. Substituting Eq. ([571) into Eqs. (IMl) and (|55|) . we 
can then fit this model to the observational data. By minimizing the corresponding total in Eq- (|17p , 
we find the best-fit parameters il-mo = 0.280, eo = —0.199 and ei = 0.214, whereas Xmin = 465.604. In 
Fig. [HI we also present the corresponding 68% and 95% confidence level contours in eo — ei plane, flmo — eo 
plane and flmo — ei plane. 



^ 0.2 
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FIG. 5: The same as in Fig. |4l except for the type II model characterized by e(a) — eo + ei{l — a). 
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Finally, we consider the logarithmic case with 



e(a) = eo + ei Ina , 



(38) 



which can be regarded as a linear expansion of e(a) with respect to the so-called e-folding time Af — Ina 
in the literature. Although e(a) diverges when a in the early universe, unlike in the same case of 
type I model, it does not cause any problem in type II model, since a"^"^*^'"^ — )• 0°° which is regular 
when a —5- 0. Substituting Eq. (|38l) into Eqs. and (PSj) . we can then fit this model to the observational 
data. By minimizing the corresponding total in Eq. dT7|) . we find the best-fit parameters fimo = 0.278, 
Co = —0.202 and ei = —0.059, whereas Xmin = 465.516. In Fig.[6l we also present the corresponding 68% 
and 95% confidence level contours in eo — ci plane, flmo — eg plane and ftmo — ci plane. 
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FIG. 6: The same as in Fig. |4l except for the type II model characterized by e(a) = eo + ei In a. 
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Best iits 


"mO ~ 0.278 


^mO = 0.281 


^^mO = 0.218 


o r* ooo 
i^mO — U.ZOO 




( \ fx oon 


iimo = 0.278 






Wj^ — — u.yoz 


"UJ^ — — U.4oo 


"UJ^ — — U.f5Df5 


eo — — U.izy 


eo — — u.iyy 








C = 2.988 


C = -3.584 


C = -2.957 


ei = 1.263 


ei = 0.214 


ei = -0.059 


Xmin 


466.317 


465.604 


563.77 


467.718 


465.635 


465.604 


465.516 


k 


1 


3 


3 


3 


3 


3 


3 


\rain / 


1.172 


1.176 


1.424 


1.181 


1.176 


1.176 


1.176 


ABIC 





11.265 


109.431 


13.379 


11.296 


11.265 


11.177 


AAIC 





3.287 


101.453 


5.401 


3.318 


3.287 


3.199 


Rank 


1 


3 ~ 4 


7 


6 


5 


3 ~ 4 


2 



TABLE I: Summarizing all the 7 models considered in this work. Here, we label the type I models characterized 
by /(a) ~ fo , fid) ~ fo + C(l ^ o,) with the condition ^ > 1 — ^^g, and /(a) ~ fo + C(l ~ o.) without the 
condition ^ > 1 — SI^q ^ IPL, ICPLw and ICPLwo, respectively. Also, we label the type II models characterized 
by e(a) = eo a'^^ , e{a) = eo + ei (1 — a), and e(a) = eo + ei In a as IIPL, IICPL and IILog, respectively. 

V. SUMMARY AND DISCUSSIONS 

In this work, we considered the cosmological constraints on the interacting dark e nerg y models. We 
generalized the models considered previously by Guo et al. [iSl, Costa and Alcaniz |l6l |. and we have 
discussed two general types of models: type I models are characterized by Pxl Pm = f{o) and /(a) can 
be any function of scale factor a, whereas type II models are characterized by p„i — Pmo a^'^+'^*^°-' and 
e(a) can be any function of a. We obtained the cosmological constraints on the type I and II models with 
power-law, CPL-like, logarithmic /(a) and e(a) by using the latest observational data. 

Some remarks are in order. Firstly, here we briefly justify the interaction forms considered in the 
present work. We take type I models as examples. For the power-law case with /(a) = fo in Eq. (|23)) . 
noting that in the case without interaction p^ oc a"^^^"*"™^-* and pm oc a~^, from definition Eq. ([9]), 
it is reasonable to parameterize /(a) = p^/p„i oc a^, where ^ measures the severity of the coincidence 
problem [13, Ell- For the CPL case with f (a) = fo + £,{1 — a) in Eq. (^5]) and the logarithmic case with 
/(a) — fo + S^lna in Eq. (P^ . noting that the Taylor series expansion of any function F{x) is given by 
F{x) = F{xo) +Fi{x- xo) + (Fa/ 2!) {x - xq)^ + (i^s/ 3!) {x - xo f + . . . , the CPL and logarithmic cases 
can be regarded as the Taylor series expansion of / with respect to the scale factor a and the e-folding 
time TV = In a up to first order (linear expansion), similar to the well-known EoS parameterizations 
w{a) = Wo + Wail — a) and w{z) = wo + wi z. 

Secondly, we would like to briefly consider the comparison of these models. For convenience, we also 
consider the well-known ACDM model in addition. Fitting ACDM model to the observational data 
considered in the present work, it is easy to find the corresponding best-fit parameter $7^0 = 0.278, 
whereas Xmm — 466.317. A conventional criterion for model comparison in the literature is Xmm/'^'^/' 
in which the degree of freedom dof ^ N — k, whereas N and k are the number of data points and the 
number of free model parameters, respectively. We present the Xmm/^o/ for all the 7 models in Table HI 
On the other hand, there are other criterions for model comparison in the literature, such as Bayesian 
Information Criterion (BIC) and Akaike Information Criterion (AIC). The BIC is defined by [33l [ssj 

BIC = -21n/:™a:i; + A:lniV, (39) 

where Cmax is the maximum likelihood. In the Gaussian cases, = —2\nCmax- So, the difference in 

BIC between two models is given by ABIC ^ Axi,„ + AfclnA^. The AIC is defined by [Ulli] 

AIC = -21n/:„„.^ + 2fc. (40) 

The difference in AIC between two models is given by AAIC ~ lS.x^nin + 2Afc. In TablelH we also present 
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the ABIC and AAIC of all the 7 models considered in this work. Notice that ACDM has been chosen 
to be the fiducial model when we calculate ABIC and AAIC. From Table [H it is easy to see that the 
rank of models is coincident in all the 3 criterions {xtunl dof ^ BIC and AIC). The ACDM model is the 
best one, whereas ICPLw model is the worst one. This result is consistent with the one obtained in 
e.g. [ll]. However, it is well known that ACDM model is plagued with the cosmological constant problem 
and the coincidence problem (see e.g. [l[). On the other hand, as mentioned in the beginning of Sec. HI 
there are some observational evidences for the interaction between dark energy and dark matter, and 
the coincidence problem can be alleviated in the interacting dark energy models. Therefore, it is still 
worthwhile to study the interacting dark energy models. 
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